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(N '. Abstract 



We perform a comprehensive analysis of the pion-photon transition form factor F n ^(Q 2 ) 
involving the transverse momentum corrections with the present CLEO experimental data, in 
■ which the contributions beyond the leading Fock state have been taken into consideration. As 

is well-known, the leading Fock-state contribution dominates of F n j(Q 2 ) at large momentum 

O 

transfer (Q 2 ) region. One should include the contributions beyond the leading Fock state in 
^p Hi small Q 2 region. In this paper, we construct a phenomenological expression to estimate the 

contributions beyond the leading Fock state based on its asymptotic behavior at Q — ► 0. Our 
present theoretical results agree well with the experimental data in the whole Q 2 region. Then, 

•rH . 

we extract some useful information of the pionic leading twist-2 distribution amplitude (DA) 
by comparing our results of F 1T1 {Q 2 ) with the CLEO data. By taking best fit, we have the DA 
moments, c^/i 2 ,) = O-OC^Iq^, a^/x 2 ,) = — 0.022^q;q 2 2 and all of higher moments, which are closed 
to the asymptotic-like behavior of the pion wavefunction. 
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I. INTRODUCTION 



The pion-photon transition form factor, which relates two photons with one lightest me- 
son, is the simplest example for the perturbative application to exclusive processes. By 
neglecting k± (the transverse momentum of the constitute quarks) relative to q± (the trans- 
verse momentum of the virtual photon) in the hard-scattering amplitude, one can obtain 



the leading Fock-state formula 



1]: 



f ^ {q2) = w 2 I^ Mx,q2 



i " 

m 



0[a s (Q 2 ), /y2 



(1) 



where Q 2 = —q 2 stands 
U = 92.4 ± 0.25 MeV 



Q 

"or the momentum transfer in the process, the pion decay constant 
and 4> n (x, Q 2 ) stands for the leading Fock-state pion distribution 
amplitude (DA) at the factorization scale fi 2 = Q 2 , which can be derived from the initial 
DA 6J x, iIq) (/io stands for some hadronic scale) through QCD evolution. The initial DA 
is defined as 



,, 2^3 f d 2 k A 



0^(x,/i o ) = -— / — —Vggfakx). (2) 

Hence, the value of Q 2 F ni {Q 2 ) tends to be a constant (2f n ) for asymptotic DA: 
(pas(x, Q 2 )|q2^ 00 = 6x(l — x). However, it has been argued that the fcj_-dependence of 
the pion wavefunction can not be safely neglected at the end-point region Xi — > 0, 1 [xi 
the momentum fraction of the constitute quarks in pion) and Q 2 ~ a few GeV 2 . Under 
the light-cone (LC) pQCD approach, the leading contribution to F nl {Q 2 ) that keeps the 
/^-corrections in both the hard-scattering amplitude and the wavefunction can be written 
as flfl: 

F^(Q 2 ) = 2V3(e 2 u - e 2 ) f\dx] f ^±S> qq (x, kj x T H (x, x', kj, (3) 

where [dx] = dxdx'5(l the quark charges in unites of e, ty qq -(x, kj_) stands for 

the leading Fock-state wavefunction and the hard-scattering amplitude T^(x,x',kj_) takes 
the form, 

= (4) 
With the help of Eqs. (13141) . Ref.[5j performed a careful analysis of the quark transverse- 
momentum effect to F 7r7 (Q 2 ). They pointed out that the transverse-momentum dependence 
in both the numerator and the denominator of the hard-scattering amplitude are of the 
same importance and should be considered consistently. Similar improved treatment has 



also been done in Refs 



11 



121 ]. It was shown that pQCD can give the correct 



prediction for the pion-photon transition form factor that is consistent with the present 
experimental data by keeping the fc^-dependence in both the hard-scattering amplitude and 
the pion wavefunction and by properly choosing of the pion wavefunction. 

It should be noted that Eqs. CTfB"]) were obtained by assuming the leading Fock-state 
dominance. This approximation is valid only for large Q 2 region and one can not expect 
that these expressions can describe the present experimental data well in low Q 2 region. 
Refs.flfl show that the approximation of the leading Fock-state dominance to the pion 
electro-magnetic form factor is valid as Q 2 ^ 4 GeV 2 , which is improved to be Q 2 ^ 1 GeV 2 
by including the next-to- leading order (NLO) contribution 15|, [l6|. A similar discussion has 
been done in Ref. 5| for the pion-photon transition form factor. These references tell us that 
one should take into account the higher Fock states' contributions as Q 2 < a few GeV 2 . 
In fact, it has been shown that the expression ([3]) gives half contribution to F nl (Q 2 ) as one 
extends it to Q 2 = [4]. It means that the leading Fock state contributes to -F^O) only half 
and the remaining half should be come from the higher Fock states. Both contributions from 
the leading Fock state and the higher Fock states are needed to get the correct tt° — > 77 
rate [4|. Any attempt that involves only the leading Fock-state contribution to explain 
both the 7T° — > 77 rate and the pion-photon transition form factor for low Q 2 region is 
incorrect. It should be pointed out that the above conclusion does not contradict with that 
of Ref. 17|, where with the help of an "effective" two-body wavefunction that includes the soft 



contributions from the higher Fock components, the authors pointed that the contributions 
corresponding to higher Fock states in a hard region appear as radiative corrections and are 
suppressed by powers of (a s /n) ~ 10%. 

In this paper, we will take the contributions from both the leading Fock state and the 
higher Fock states into consideration. Especially, we will discuss how to consider the con- 
tributions beyond the leading Fock state at low Q 2 region and give a careful analysis of the 
pion-photon transition form factor in the whole Q 2 region. Furthermore, we can learn more 
information of the leading twist-2 DA from the present CLEO data, since the pion-photon 
transition form factor in the simplest exclusive process only involves one pion and the con- 
tributions from the higher twist structures and higher helicity states are highly suppressed 
(at least by 1/Q 4 ) in comparison to the leading twist-2 pion wavefunction. In the present 
paper, we will not include the NLO contribution into our formulae. Since we need a full 
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NLO result in order to be consistent with our present calculation technique, in which the 
effects caused by the transverse-momentum dependence in the hard-scattering amplitude 
and the pion wavefunction and by the Sudakov factor should be fully considered, however 
such a full NLO calculation is not available at the present 1 . 

The paper is organized as follows. In Sec. II, we analyze the contributions to -F 7r7 (Q 2 ) 
beyond the leading Fock state at low Q 2 region under the LC pQCD approach and give 
a complete expression for F n .y(Q 2 ) in the whole Q 2 region. In Sec. Ill, we discuss what 
we can learn of the pionic leading Fock-state wavefunction/DA in comparison with CLEO 
experimental data. Some further discussion and comments are made in Sec. IV. The last 
section is reserved for a summary. 



II. AN EXPRESSION OF F^{Q 2 ) FROM ZERO TO LARGE Q 2 REGION 



First, we give a brief review of the LC formalism l|, [19|, |20|]. The LC formalism provides 
a convenient framework for the relativistic description of hadrons in terms of quark and 
gluon degrees of freedom and for the application of pQCD to exclusive processes. The LC 
Fock-state expansion of wavefunction provides a precise definition of the parton model and 
a general method to calculate the hadronic matrix element. As for the pion wavefunction, 
its Fock-state expansion is 



k) = E + E l?&>*«© + ■ 



(5) 



where the Fock-state wavefunctions ^ n (xi, k_u, A») (n = 2, 3, • • •) satisfy the normalization 
condition 



E J[dxd 2 k ± } n \^ n (xi,k A 

n,\i 



i, Aj) 



16vr 3 



(6) 



. Aj is the helicity of the 



with [dxd 2 k ± } n = 16vr 3 ,5(l-Er=i^)^(E? = ik ±J )nr=i 
constituents and n stands for all Fock states, e.g. ^2 — ^qq- It should be pointed out that 
we have J[dxd 2 k±} 2 \^qq(xi, k_u, Aj)| 2 < 1 for the leading Fock state. 

The pion-photon transition form factor is connected with the 7r°77* vertex in the ampli- 



1 The present NLO result is derived without considering the transverse-momentum dependence in the hard- 
scattering amplitude and the wavefunction, see e.g. 
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FIG. 1: Typical diagrams that contribute to the pion-photon transition form factor i ? 7r7 (Q 2 ), where 
x' = (1 — x). The rightmost shaded oval with a slant pattern stands for the strong interactions. 

tude of en — > as 

T M = -ie 2 F^{Q 2 )e^ aP P»e a q^ (7) 

where P and q are the momenta of the incident pion and the virtual photon respectively, 
and e is the polarization vector of the final on-shell photon. To simplify the hard-scattering 



amplitude, we adopt the Drell-Yan-West assignment at the 'infinite-momentum' frame [211 ] : 
q = (g+,<Z-,q±) = (0,|qJ 2 /P+,q ± ) and P = (P+,P~,P ± ) = (1,0, ± ), where P+ is 
arbitrary because of Lorentz invariance and it is taken to be 1 for convenience. Then F nl is 
given by [l| 

F^{Q 2 ) = ./ + v v (8) 

where Q 2 = -q 2 = \q±\ 2 , e = (0, 0, e±), e± ■ q ± = and e± x q ± = e_Lig_L 2 - e_L2<?_Li- 

As illustrated in Fig.([Tj), there are two basic types of contribution to F ni (Q 2 ) , i.e. 
F^\Q 2 ) and F^ V \Q 2 ). Fff(Q 2 ) comes from Fig.fllk), which involves the direct anni- 
hilation of (gg)-pair into two photons, i.e. the leading Fock-state contribution that dom- 
inates the large Q 2 contribution. F^ V \Q 2 ) comes from Fig.fPJ)), in which one photon 
coupling 'inside' the pion wavefunction, i.e. strong interactions occur between the photon 
interactions that is related to the higher Fock states' contributions. An interpretation for 
F^ V \Q 2 ) can be given under the operator product expansion (OPE) approach [171 ] . Under 
the OPE approach [221], the nonperturbative aspects of the hadron dynamics are described 
by matrix elements of local operators. In particular, the longitudinal momentum distribu- 
tion is related to the lowest-twist composite operators. While by taking into account the 
transverse-momentum effects, one needs to consider matrix elements of higher-twist com- 
posite operators in which some of the covariant derivatives appear in a contracted form like 



D 2 = D^D^. By using the equation of motion for the light quark, ^D^q = 0, one can 
convert a two-body quark- ant iquark operator g{7 A11 D M2 . . .D pkn }D 2 q into the "three-body" 
operator qij^D^ . . . D^ n }{a^ u G^ u )q with an extra gluonic field G^ v being involved, which 
is just related to the higher Fock state of pion. 

The first type of contribution F^'(Q 2 ) stands for the conventional leading Fock-state 
contribution. Under the LC pQCD approach and by keeping the full fey-dependence in 
both the hard-scattering amplitude and the wavefunction, the expression for F^(Q 2 ) is 
the one that is given in Eqs. (131141) . where the terms involving the higher helicity states 
(Ai + A2 = ±1, Aj stands for the corresponding helicity of the two constitute quarks of pion) 
and the higher twist structures of the pion wavefunction are not explicitly written. Since 
by direct calculating, one may observe that the contributions from the higher helicity states 
and the higher twist structures of the pion wavefunction are suppressed by at least 1/Q 4 to 
that of the usual helicity state (Ai + A2 = 0) of the leading Fock state, which agrees with 
the discussion made in Ref.[9| 2 . 

As for the second type of contribution Ff£j v \Q 2 ) , it is difficult to be calculated in any Q 2 
region. If treating the photon vertex in Fig. ([lb) as a vector meson dressed photon vertex, 
Fig. (Jib) can be calculated approximately under the vector meson dominance (VMD) ap- 



proach, see e.g. Ref. }23| for a review and Refs. 24], |25( for an explicit VMD calculation of the 



pion electro-magnetic form factor. By adopting the VMD approach to approximate Fig. ([lb), 
one needs to introduce some undetermined coupling factor either for VMD1 or VMD2 for- 



mulation [23||, which together with the undetermined parameters in the pion wavefunction 
can not be definitely determined by the CLEO experimental data of the pion transition 
form factor only and some other constraints should also be taken into consideration, e.g. 
the constraint from the experimental value for the pion charge radius or the constraint from 
the experimental value for the pion electro- magnetic form factor. In fact, one usually takes 
the value of i^ 7 (Q 2 ) = F^(Q 2 ) + F^ V \Q 2 ) derived from the VMD approach to be in a 



simple monopole form |26|], i.e. F WJ (Q ) = 1/ An f n (l +Q /m 2 ) , with the p-meson mass 
m p serves as a parameter determined by the pion charge radius. For the purpose of extract- 
ing some useful information of the pion wavefunction from the CLEO experimental data, 



2 The present condition is quite different from the case of pion electro-magnetic form factor, where the 
contributions from the higher helicity states and higher twist structures are only suppressed by 1/Q 2 and 



then they can provide sizable contribution in the intermediate Q 2 region 37J, |38| . 
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we adopt the method raised by Ref.J^j to deal with F.^ V \Q 2 ) 3 . 



As stated in Ref.[4J, around the region of Q ~ 0, since the wavelength of the photon 
'inside' the pion wavefunction ~ is assumed to be much larger than the pion radius 

1/A (A is some typical hadronic scale ~ 1 GeV), we can treat such photon (nearly on-shell) 
as an external field which is approximately constant throughout the pion volume. And 
then, a fermion in a constant external field is modified only by a phase, i.e. Sa{x — y) — 
e -ie(y-x)-Ag F ^ x _ ^ Consequently, the lowest qq- wavefunction for the pion is modified 
only by a phase e~ iey ' A , where y is the gg-separation. Transforming such phase into the 
momentum space and applying it to the wavefunction, the second contribution Ff^ v \Q 2 ) 
at qj_ — > can be simplified to 



2 



v/3Q2 J 1 J J 16tt 3 [ (x'q ± 







+ (x <-> x'\ 



' (9) 

where [dx] = dxdx'8(l — x — x') and k± = |kj_|. Eq.Q gives the expression for F^ V '(Q 2 ) 
at Q 2 — > 0. Here different from Ref.jj], all qj_-terms that are necessary to obtain its first 
derivative over Q 2 are retained and the relation (e_i_ x q_i_)(kj_ x qj_) = Q 2 (e±-k±) is implicitly 
adopted. After doing the integration over k^, one can easily find that 

F$ V \0) = F%\0) = ^-L- 2 J dx%,(x, ± ), (10) 

which means that the leading Fock state contributes to i^ 7 (0) = F^(0) + F^ v \0) only 
half, and one can get the correct rate of the process 7r° — ► 77 provided that the two basic 
contributions F^p(0) and F^ v \0) are considered simultaneously. By taking into account 
the PCAC prediction Q, F^ 7 (0) = l/(4vr 2 ^), one can obtain the important constraint of 
the pion wavefunction, i.e. 

.' 1 dx* w -(x,k J . = 0) = ^. (11) 

Without loss of generality, we can assume that the pion wavefunction depending on k^ 
through k\ only, i.e. ^ q g(x,k ± ) = ^ q q(x,k 2 ± ) 4 . Then Fff(Q 2 ) (Eq.©) can be simplified 



3 A careful VMD calculation of the pion transition form factor along the line of Refs.(24|, [25| can be used 
as a cross check of our results, however it is out of the range of the present paper. 
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The spin-space Wigner rotation might change this property for the higher helicity components as shown 
in Ref. 28|. Since the higher helicity components' contribution are highly suppressed for the present case, 



we do not take this point into consideration in the present paper. 
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after doing the integration over the azimuth angle as 9] 

1 r 1 dx 



V qq (x, ki)dk\ 



As/?>ii 2 Jo xQ 2 Jo 
Similarly, for the first derivative of F^ V \Q 2 ) over Q 2 , we have 



fS v) \Q 2 )V 



1 



8^ 



7l z 



a 



1 rx 



'Vqg(x, hi) 



docdk i 



(12) 



(13) 



dQ 2 Jo Jo V x2 Q 2 
Furthermore, the leading twist-2 pion DA at the factorization scale \x can be simplified as 

y/3 r» 



8tt 2 U Jo 



With the help of Eqs.(|12][14j), Fff(Q 2 ) can be rewritten as 

\2\ dx 



(14) 



(15) 



T 7 ^ J 3 Q 2 

Note that Eq.( fl5i) is different from Eq.(TjQ) only by replacing <p n (x,Q 2 ) to 4> n (x,x 2 Q 2 ). It 
means that the leading contribution to F nl (Q 2 ) as shown in Eq.([3]), which was given by 
keeping the /c^-corrections in both the hard-scattering amplitude and the pion wavefunction, 
can be equivalently obtained by setting the upper limit for the integral of the pion DA to 
be [/z 2 = x 2 Q 2 ]. The x-dependent upper limit [x 2 Q 2 ] affects F^\Q 2 ) from the small to 
intermediate Q 2 region, and such effect will be more explicit for a wider pion DA, such as 
the CZ (Chernyak-Zhitnitsky)-like model [s| that emphasizes the end-point region in a strong 
way, as has been discussed in Ref.^. In the literature, the pion DA is usually expanded in 
Gegenbauer polynomial expansion as 

oo 

i + E^(^ 2 )^ 2 (0 



(/>„(x,tl ) = <fia S (x) ■ 



n=l 



(16) 



where £ = (2x — 1), C^/ 2 (£) are Gegenbauer polynomials and a 2n (/i 2 ), the so called Gegen- 
bauer moments, are hadronic parameters that depend on the factorization scale fi. The 
Gegenbauer moments a2„(/i 2 ) can be related to a^n^o) with the help of QCD evolution, 
where Up stands for some fixed low energy scale. To leading logarithmic accuracy, we 



have 



10 



a 



(17) 



where (3q = 11 — 2nj/3, a s (Q 2 ) = 47r/[/3 \n(Q 2 / Aq CD )} and the one-loop anomalous dimen- 
sion is 

%C F L{2n + 2)+ lE --- - — *- — ). (18) 



l { o 2n) 



4 (2n+l)(2n + 2), 
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We need to know 4> 1T (x,x 2 Q 2 ) and F^ V \Q 2 ) to get the whole behavior of F 7ry (Q 2 ), i.e. 

^(Q 2 )=Fi v 7 \Q 2 ) + F^ v \Q 2 ), (19) 

where F^(Q 2 ) is determined by (j) n (x, x 2 Q 2 ). The (j) w (x, x 2 Q 2 ) depends on the behavior 
of ty q q(x, kj_) and its Gegenbauer moments can not be directly obtained from the QCD 
evolution equation (1TT1) . since [x 2 Q 2 ] can be very small and then the Landau ghost singularity 
in the running coupling a s can not be avoided. As for F^ V \Q 2 ), Eq.(|9]) presents an 
expression only at Q 2 ~ region, and it can not be directly extended to the whole Q 2 
region. Ref. 9| did an attempt to understand the higher Q 2 behavior of F 1X1 {Q 2 ) within 
the QCD sum rule approach, i.e. they raised a simple picture: the sum over the soft 
qG ■ ■ -Gq Fock components is dual to gg-state generated by the local axial vector current. 
Furthermore, they raised an 'effective' two-body pion wavefunction that includes all soft 
contributions from the higher Fock states based on the QCD sum rule analysis and then 
calculated F nl (Q 2 ) within the pQCD approach. Here we will not adopt such an 'effective' 
pion wavefunction to do the calculation, since we plan to extract some information for the 
leading Fock-state wavefunction by comparing with the CLEO experimental data. 

In order to construct an expression of F^ V \Q 2 ) in the whole Q 2 region, we require the 
following conditions at least: 

i) F^ V \Q 2 )\ Q2=0 should be given by Eq.fllOj). 

ii) F^)'(g 2 )| Q2 ^ = dF^ v \Q 2 )/dQ 2 \ Q 2^ should be derived from Eq.©. 

iii) ^j^l -> 0, as Q 2 - oo. 

One can construct a phenomenological model for F^ V \Q 2 ) that satisfies the above three 
requirements. It is natural to assume the following form 

F^ V \Q 2 ) = (1 + g2//?2)2 , (20) 
where k and a are two parameters that can be determined by the above conditions (i, ii) , i.e. 

a = -F^JO) = (21) 

and 



jjmO) /no's 

As for the phenomenological formula (J2"0"j) . it is easy to find that F^ V \Q 2 ) will be suppressed 
by 1/Q 2 to F^XQ 2 ) in the limit Q 2 -> oo. Such a 1/Q 2 -suppression is reasonable, since the 
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FIG. 2: The fitting curve (the solid line) for Q 2 F. K1 (Q 2 ) from the CLEO and CELLO experimental 



data 



30, 



311 ] . where a shaded band shows its uncertainty of ±20%. 



phenomenological expression (1201 can be regarded as a summed up effect of all the high twist 
structures of the pion wavefunction, even though each higher twist structure is suppressed 
by at least 1/Q 4 |9j. 



III. CALCULATED RESULTS WITH THE MODEL WAVEFUNCTION 



The CLEO collaboration has measured the 77* — > 7r° form factor 30] . In this experiment, 
one of the photons is nearly on-shell and the other one is highly off-shell, with a virtuality 
in the range 1.5 GeV 2 - 9.2 GeV 2 {jjoj]. There also exists older experimental results obtained 
by the CELLO collaboration 3jj . By comparing the theoretical prediction with the experi- 
mental results, it provides us a chance to determine a precise form for the leading Fock-state 
)ion wavefunction. Similar attempt to determine the pion DA has been done in literature 



32 



33| . e.g. Ref.[32j used the QCD light-cone sum rule analysis of the CLEO data to obtain 



parameters of the pion DA. In FigfSJ we show the fitting curve for Q 2 F nl (Q 2 ) (derived by 
using the conventional x 2_n tting method described in Ref. 34j with slight change to make 
the curve more smooth) from the CLEO and CELLO experimental data, i.e. under the 
region of Q 2 G [0.5, 10.0] GeV 2 , Q 2 F nl (Q 2 ) ~ [8.81 x 10- 7 (Q' 2 ) 5 - 4.78 x 10- 5 (Q' 2 ) 4 + 9.96 x 
10~ 4 (Q' 2 ) 3 - 1.01 x 10~ 2 (Q' 2 ) 2 + 5.29 x 10~ 2 (Q' 2 ) + 4.48 x 10~ 2 ] GeV with the dimensionless 
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parameter Q" 2 = Q 2 /GeV 2 . Here the shaded band shows its ±20% uncertainty 5 . In fact, 

nnnnnnnri 

most of the results given in literature, e.g. Refs. [5|, la, |7|, l8|, |9|, [1Q|, llll, |32|, are mainly within 



such region. The shaded band (region) for Q 2 F nl (Q 2 ) can be regarded as a constraint to 
determine the pion wavefunction, i.e. the values of the parameters in the pion wavefunction 
should make Q 2 F nl (Q 2 ) within the region of the shaded band as shown in Fig.(T5]). 

Now we are in position to calculate the pion-photon transition from factor with the help of 
Eq. ffT9l . As has been discussed in the last section, we need to know the leading Fock-state 
pion wavefunction ^ g q(x, kj_) so as to derive ^(x, x 2 Q 2 ) that is necessary for F^(Q 2 ) 
and to derive the values of a and k for F^ V \Q 2 ). Several non-perturbative approaches 
have been developed to provide the theoretical predictions for the hadronic wavefunction. 
One useful way is to use the approximate bound-state solution of a hadron in terms of the 
quark model as the starting point for medeling the hadronic wavefunction. The Brodsky- 
Huang-Lepage (BHL) prescription [4j] for the hadronic wavefunction is obtained in this way 
by connecting the equal-time wavefunction in the rest frame and the wavefunction in the 
infinite momentum frame. In the present paper, we shall adopt the revised LC harmonic 
oscillator model as suggested in Ref. 28j to do our calculation, which is constructed based on 
the BHL-prescription. As discussed in the above section, the contribution from the higher 
helicity states (Ai + A 2 = ±1) is highly suppressed in comparison to that of the usual helicity 
state (Ai + A2 = 0), so we only write down the form of the pion wavefunction for the usual 
helicity state: 

Vqq(x,k±) = <Pbhl(x, k ± )x K (x,k ± ) = Aexp " g^^ _ ~\ X (x,^±), (23) 
with the normalization constant A, the harmonic scale (3 and the quark mass m to be 



determined. The spin-space wavefunction % (x, kj_) can be written as |28|, x K (x,k 



m/^Jm 2 + k\ with k± = |kj. By taking the BHL-like wavefunction ([23]), F^(Q 2 ) (Eq.(fT2l)) 
can be simplified as 



1 , Am/3 

(IX < — = ; 

I VQk 3/2 Q 2 



Erf 



y/m 2 + x 2 Q 2 
2(3^ 



xx 



Erf 



77T 



2(3 y/2 



xx 



(24) 



where the error function Erf (a;) is defined as Erf (a;) = Jq e f2 dt. And similarly, for the 
limiting behaviors of F^ V \Q 2 ) that are necessary to determine the parameters a and n, 



5 It is so chosen since the sum of the statistical and systematic errors of the experimental data is < 
±20% Si- 
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FIG. 3: The curve for the value of m versus (3, which shows that if m is in the reasonable region 
[0.20 GeV,0.40 GeV] then (3 G [0.48 GeV,0.70 GeV]. 



we obtain 



A 



8V3tt 2 Jo 



exp 



m 



1 >P 2 xx' 



dx 



and 



-A 



x 

/ — (m 2 + Axx' ' p 2 ) exp 
Jo X 



m 



i(3 2 xx' 



(25) 



dx, (26) 



where x' — 1 — x. 

The possible range for the parameters of the pion wavefunction (1231) can be derived by 
comparing our result of Q 2 F nl (Q 2 ) with the experimental data. Further more, the model 
wavefunction should satisfy the following two conventional constraints: 

• The pion wavefunction satisfies the normalization that is derived from the tt — > fiu 



process 



2Sj: 



dx 



d 2 k. 



^na(x, k| 



(27) 



io J\* x \* <(t > 16^ ^'""^ 2^3' 
where the cut |kj_| 2 < fi 2 has been explicitly included. Substituting the pion wave- 
function fl23|) into Eq. (l27"|) . we obtain 



i Amj3^x{l - x) 
4\/27r 3 / 2 



Erf 



m 2 + \x 2 
\ 8/? 2 x(l - x) 



Erf 



\ 8/? 2 x(l 



<ix 



2y/3' 



(28) 

As for the model wavefunction fl23l) . it can be found that the contribution from higher 
|kj_ | region to the wavefunction normalization drops down exponentially, e.g. by taking 
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reasonable values for the wavefunction parameters, the wavefunction in the region of 
| kj_ | > 1 GeV only contributes about 5% to its total normalization, which changes to 
be less than 1% in the region of |kj_| > 2 GeV. For clarity, we take the factorization 
scale /i to be /i = /io — 2 GeV, where /xo stands for some hadronic scale that is of order 
(9(1 GeV) and the choice of /xq — 2 GeV is also close to the virtuality of the photons 
in the central region of the experimental data. By taking 1 GeV < xxo < 2 GeV that 
is of order /xo ~ 0(1 GeV), we find that the following results will be slightly changed, 
especially for the DA moments due to the fact that they satisfy the QCD evolution 
equation (1T7I) within errors. A similar discussion on this point can also be found in 



Ref. 



32J. Furthermore, one can safely take /x = oo to simplify the calculation due to 



the smallness of the wavefunction in the region of [|kj_| > /Xq = 2 GeV]. 



Another constraint, as shown in Eq. ffTTj) . for the pion wavefunction can be derived 
from 7T° — > 77 decay amplitude [4], which can be further simplified as 



o 



l 

Aexp 



2 

m 



dx = ^. (29) 



8(3 2 x(l-x 

Solving Eqs.( |28|29l) numerically, we obtain an approximate relation for m and (3, i.e. 

6.00^ = 1.12 + 1.3l) (j + 5.47 x 10 1 ) , (30) 

which shows that the value of m is decreased with the increment of (3. More explicitly, 
Fig. ([3]) shows that (3 should be within the region of [0.48 GeV, 0.70 GeV] so as to restrict m 
within the reasonable region of [0.20 GeV, 0.40 GeV]. 

Before doing the numerical calculation for F^IQ 2 ), we note a naive interpolation formula 
for both perturbative and non-perturbative regions has been proposed by Brodsky and 



Lepage [lj, which is similar to the monopole form derived from VMD approach 26], i.e. 



F - m = «S) t 1 - ^ ) • s ° = = 67 GeV " ~ m > (31) 



where the NLO perturbative contribution (_| Qs (Q ) \ been added to the original result 

nn 

according to the suggestion of Ref. [9|, H8j . 

The conventional value for the constitute quark mass m is around 0.30 GeV. In the 
present section, we concentrate our attention on the case of m — 0.30 GeV and we will 
study the uncertainty caused by varying m within a wider region of [0.20 GeV, 0.40 GeV] 
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Q 2 (GeV 2N 



FIG. 4: Q 2 F 7T ~ f (Q 2 ) that is defined in Eq. ()19p and is calculated with the model wavefunction (|23|) 
for the case of m = 0.30 GeV ((3 = 0.55 GeV), which is shown by a dash-dotted line. The dashed 



line is best fit for the CLEO experimental data 



30, 



31 ], the solid line is from the interpolation 



formula (J3TJ) and the dotted line is the results for Q 2 F^\Q 2 ) only. 



in the next section. We show the pion-photon transition form factor Q 2 F nl (Q 2 ) with the 
model wavefunction (1231) for the case of m = 0.30 GeV in Fig.®, where the values for 
P = 0.55 GeV and A = 25.6 GeV^ 1 can be obtained by using Eqs. fl271l30|) . Fig.® shows 
that by taking the model wavefunction (|23l) with m = 0.30 GeV, the predicted value of 
Q 2 F nl (Q 2 ) agrees well with the experimental data. We also show the leading Fock-state 
contribution Q 2 F^\Q 2 ) in Fig.®, which is drawn in a dotted line and its value is lower 
than the experimental data especially in low Q 2 region. As a comparison, we draw the curve 
for the interpolation formula fl3T|) in Fig.®, where similar to Ref.[3j|, we adopt the one 
loop a s -running with A.qq^ = 312 MeV to do our calculation. It shows that Q 2 F^(Q 2 ) 
agrees with the data especially for the low Q 2 region, which is reasonable as the effective 
value of so in F^ L (Q 2 ) is determined by the known behavior of Q 2 — > 0. 

IV. DISCUSSION AND COMMENT 

A. Information of the leading Fock state 



As shown in Fig.®, F^^Q 2 ) agrees well with the experimental data by taking the leading 
Fock-state pion wavefunction (|23|) with m = 0.30 GeV. We also show the leading Fock-state 
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contribution Q 2 Ff^\Q 2 ) in Fig.(jlj), which is lower than the experimental data in low Q 2 
region. This shows that one should take into account the higher Fock states in small to 
intermediate Q 2 region. In fact, it has been found that the leading Fock-state contribution 
F^(Q 2 ) fail to reproduce the Q 2 = value corresponding to the axial anomaly 5|, |7|, i.e. 
it gives only half of what is needed to get the correct 7r° — > 77 rate [271]. And to make 
a compensation, in Refs. , the leading Fock-state contribution F^{Q 2 ) has been 

enhanced by replacing the leading Fock-state wavefunction to an 'effective' valence quark 
wavefunction that is normalized to one. By taking the 'effective' pion wavefunction with 
the asymptotic-like DAs, the authors found an agreement with the experimental data for 
F nl (Q 2 ) [6, [9I, 17]. However, such an 'effective' pion wavefunction is no longer the leading 



Fock-state wavefunction itself and the probability of finding the leading Fock state in pion 
should be less than one. 

By substituting the pion wavefunction ( |23i) into the pion electro-magnetic form factor, 
one can obtain some useful information, such as the probability of finding the leading Fock 



state in pion P qq , the mean square transverse-momentum of the leading Fock state (k 2 ) 



qq 



and the charged mean-square-radius (r 2 + ) qq . In Refs. 37J, |39j, the authors have done such 
a calculation within the LC pQCD approach. By adopting the formulae derived in Ref. 37 ] 
(Eqs. (24,25,27) there), one may obtain: 



P qq = 56%, (kl) M - 



(0.502 GeV) 2 , (r 2 +) 9 « = (0.33 fm) s 



(32) 



where we have taken m = 0.30 GeV. All the values in Eq. (1321) are summed results for all 
the helicity states Ai + A2 = (0, ±1) of the leading Fock state. Eq. (1321) shows that the 



value of (r 2 +) qq is smaller than the pion charged radius (r 2 )1 xpt = (0.671 ± 0.008 fm) 2 (40], 



but it is close to the value as suggested in Refs.|37j, |4jJ. Such small (r 2 +) qq for the leading 
Fock-state wavefunction is reasonable, since the probability of leading Fock state P q q is only 
56%, which confirms the necessity of taking the higher Fock-states into consideration to give 
full estimation of the pion electromagnetic form f actor /pion-photon transition form factor, 
especially for lower Q 2 regions. 
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FIG. 5: Q 2 F W1 (Q 2 ) with the phenomenological model (JH for F^ V \q 2 ), with the upper edge 
of the band for m = 0.40 GeV and the lower edge of the band for m = 0.20 GeV. The dashed 
line is the best fit of Q 2 F 7r -y(Q 2 ) for wavefunction (|23p with m = 0.30 GeV. As a comparison, the 
interpolation formula (|3ip is shown by a solid line. 

B. Moments of the DA from the CLEO data 



In this subsection, we take a wider region for m, i.e. m = 0.30loio GeV, to study the 
Gegenbauer moments of pion DA. Under such region for m, we show the value of Q 2 F 7Ty (Q 2 ) 
in Fig. (pj]). The value of Q 2 F nl (Q 2 ) will increase with the increment of m and the uncertainty 
caused by varying m within the region of [0.20 GeV, 0.40 GeV] is about ±20%. By comparing 
Fig.flHJ) with Fig.([2J), it can be found that m G (0.20 GeV, 0.40 GeV) is a reasonable region 
for Q 2 F 7TJ (Q 2 ), which is in agreement with the experimental data in the whole Q 2 region. 

The pion DA in Eq.(j2J) can be defined through the matrix elements of local operators 
between the physical pion and the vacuum state 
choosing the frame in which P = (1, 0, OjJ that has been adopted in Sec. II to calculate the 
pion-photon transition form factor, its explicit form is given by the following equation 441 ] : 

\ (Q 2 ) 



42 



431 ] . In the light-cone framework, by 



dz 
~2T 



D i(x\-X2)z /2 



x ( 



■it 



2v/2 



-z, z\d(z) 



tt(P) 



(33) 



=z±=0 



where the factorization scale is taken to be Q 2 and the path-order factor [—z, z] = 
V exp (ig s j\ dsA + (zs)z~ /2^j stands for the conventional Wilson line connecting the points 
—z and z. The line integral vanishes in the light-cone gauge due to A + = 0. The behavior 
of <fr(xi, Q 2 ) for the large Q 2 is dominated by the behavior of the operator T(u(—z)d(z)) for 
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z 2 = 0(1/ Q 2 ). One can apply the usual OPE and only gauge invariant operator that is 
actually contributes to the matrix element in the light-cone gauge. Therefore only the qq 
component is required in the above definition for the leading power behavior 43j. Also the 
corresponding higher twist operators of the wave function matrix elements are suppressed 
at the short distance by power of l/Q 2 - Thus, the DA (p(xi,Q 2 ) defined in Eq. (l3"3l is the 
amplitude for finding the qq component that is collinear up to the scale Q 2 . On the other 
hand, the matrix element (0|'U7 Ai 75(i|7r(P)) can be expressed by the following expansion 45] 



(0\u(-z)^ 5 [-z,z}d(z)\n(P))^ = 

iV2UP, f 1 due^ 2u -V p - z M%Q 2 ) + A=rnl-£- t due^-^g^Q 2 )^) 
Jo v2 -r • z Jo 

where 4> n is the leading twist-2 DA, g„ is the twist-4 DA. Under the light-cone gauge, the 
DA 4>(xi, Q 2 ) defined in Eq. (l33j) corresponds to the leading twist-2 DA <p n defined in Eq.(!34j). 
except for an overall normalization factor. Both definitions are consistent with each other. 

We can derive the leading twist-2 pion DA that is in the usual helicity (Ai + A 2 = 0) from 
Eq.([23D, i.e. 



Am(3V3Jx(l - x) 



2^2/^3/2 



^Erf 




m 2 + /ig 


-Erf 




m 2 


) 


A 


8p 2 x(l - x) 


.\ 


8p 2 x(l-x) 





(35) 



where the fixed low energy scale /io is taken to be 2 GeV and 4> n (x, nl) satisfies the nor- 
malization Jq 1 (p^x, fil)dx = 1. By expanding Eq. fl35|) into the Gegenbauer polynomials, for 
m E [0.20 GeV,0.40GeV], we obtain 



-0 022 +a026 



a 2 (^) = 0.002™, aM) 
aM) = -0.014^, a 8 (/$ = -0.006™, 



(36) 



where the center value is for taking m ~ 0.30 GeV that best fits the CLEO experimental data, 
i.e. it has the minimum x 2 -value, and the ellipsis stands for higher Gegenbauer moments. For 
the values of a2 n (/x 2 ) in other factorization scales, they can be derived by QCD evolution e.g. 
Eq. ffTTl) . Eq. fl3"6"|) shows: A) the leading twist-2 pion wavefunction (12"3"|) is asymptotic-like, 
since a2 n (n > 1) are much smaller than ao = 1. More explicitly, the first inverse moment of 
the pion DA at energy scale /zo, Jo dx(j) 7T (x, nl) /x = 3(1 + a 2 + 04 + + a$) e (2.71,3.20), 
which is near the same value as for the asymptotic wavefunction with a2 n = (n > 1). Such 
a conclusion for pion wavefunction agrees with that of Ref. 46j and also agrees with a recent 
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study that is based on the nonlocal chiral-quark model from the instanton vacuum 47j . B) 
02) fl 4 w ih increase with the decrement of m, and 02 > if m < 0.30 GeV. C) the absolute 
values of a 4 , a 6 and a 8 are comparable to a 2 for bigger m (e.g. m ~ 0.30 GeV); but they 



are suppressed to 02 about one order for smaller m 



have been studied in various processes, cf. Refs.[15l. 
The lattice result of Ref. 



'he value o: 



32 



35 



36 



the Gegenbauer moments 



47, 



48 



49 



50 



51 



52 



53]. 



50 



521 ] prefers a narrower DA with 02(1 GeV 2 ) = 0.07(1), while the 
5l| prefer wider DA, i.e. they obtain a 2 (l GeV 2 ) = 0.38 ± 0.23l(^ and 



lattice results 

a 2 (l GeV 2 ) = 0.364 ±0.126 respectively. However as argued in Ref. 



54j , the accuracy of the 



lattice results needs to be further improved, e.g. the results in Ref. [511] for the second moment 
(£ 2 ) are obtained for the "pion" with the masses fi^ > 550 MeV and then extrapolated 



to the chiral limit fi n — > 0. These re 
the most recent one is done by Ref. 



erences favor a positive value for 02(1 GeV 2 ) and 
49j, which shows that a 2 (l GeV 2 ) = 0.19 ± 0.19 and 
04(1 GeV 2 ) > —0.07 by analyzing the leptonic mass spectrum of B — > ixlv. Here, the range 
of m should be reduced tome [0.20 GeV, 0.30 GeV] if we require 02(1 GeV 2 ) > for the 
pion DA. Or inversely, we have 02(1 GeV 2 ) £ [0,0.08] with the help of the QCD evolution 
equation (IT?]) for m £ [0.20 GeV, 0.30 GeV]. 



C. Comparison with the broad wavefunction 

One typical broad wavefunction is described by the CZ-like wavefunction, which can not 
be excluded by the pion-photon transition form factor [s| although it is disfavored by the 



pion structure function at x — > 1 



28] 



We take the CZ-like wavefunction as 



Vgg Z (x, kj.) = A{1 - 2xf exp 



k 2 | + m 2 



X K (x,kj 



(37) 



i(3 2 x(l -x)_ 

where an extra factor (1 — 2x) 2 is introduced into the pion wavefunction ( 1231) jjj. Following 
a similar procedure, we find that the relation between m and (3 changes to 



6.00^1.24(1 + 2.12) (I+4.58X10 1 ). 



(38) 



Similarly, for m = 0.30 GeV, we have 



P qq = 73%, (ki) M - = (0.496 GeV)^, (r 2 + ) w " = (0.45 fm) 2 , 



(39) 
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FIG. 6: Comparison of Q 2 F 7T y(Q 2 ) that is derived from two different types of wavefunctions, i.e. 
the BHL-like wavefunction (the solid line) and the CZ-like wavefunction (the dash-dotted line), 
under the condition of m = 0.30 GeV. 

where all the values are summed results for all the helicity states Ai + A 2 = (0, ±1) of 
the leading Fock state. One may observe that the second Gegenbauer moment is always 
dominant over other higher Gegenbauer moments for the CZ-like DA. And for the first inverse 
moment of the CZ-like pion DA at energy scale /io, we obtain f dxcj)^ z (x, fjfy/x = 4.69. 

In Fig. (|6]), we make a comparison of Q 2 F nl (Q 2 ) that is derived from two different types 
of wavefunctions, i.e. the BHL-like wavefunction and the CZ-like wavefunction, under the 
same value of m — 0.30 GeV. Both the BHL-like wavefunction and the CZ-like wavefunction 
lead to Q 2 F ni (Q 2 ) within the possible region of the experimental data as shown in Fig.(f2]). 
However, the value of Q 2 F 1T ~ i (Q 2 ) derived from the BHL-like wavefunction is better than 
that of the CZ-like wavefunction. One may observe that the value of Q 2 F nl (Q 2 ) caused by 
the CZ-like wavefunction shall increase with the increment of m, so the CZ-like model can 
give a better result for Q 2 F nl (Q 2 ) in comparison to the experimental data only by taking a 
bigger value for m, e.g. at least, m = 0.40 GeV. 

As a summary, the main differences for the BHL-like wavefunction and the CZ-like wave- 
function are listed in the following: 

• By comparing with the experimental data for Q 2 F ni (Q 2 ), one may observe that m = 
0.30j^o'io GeV is a reasonable region for the BHL-like wavefunction, where the best fit 
to the experimental data is achieved when m ~ 0.30 GeV; while for the case of the 
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FIG. 7: The value of R = —rrr, — — — -mrn versus Q for the BHL-like wavefunction, whose 

value increases with the increment of m. The solid, dashed and dotted lines are for m = 0.20GeV , 
m = 0.30GeU and m = OAOGeV, respectively. 

CZ-like wavefunction, such region is shifted to a higher one, i.e. m = 0.40loio GeV, 
and the best fit to the experimental data is achieved as m w 0.40 GeV, which is 
somewhat bigger than the conventional value for the constitute quark mass of pion. 

• The difference among the Gegenbauer moments e.g. a 2 (/ig) are big due to the different 
behavior of the two models. Under the condition of m = 0.30 GeV, the two Gegenbauer 
moments a 2 (/•*<)) = °- 002 and a ^o) = -0.022 for the BHL-like model ([20]); while for 
the CZ-like model fl37D, a 2 {fj, 2 ) = 0.678 and a^/ig) = -0.024. 

• The first inverse moments are different. Under the condition of m = 0.30 GeV, for the 
case of the BHL-like model ([20]), dx<j)^ HL (x, nl)/x = 2.88, which is close to that of 
the asymptotic DA; while for the case of the CZ-like model (1371) . Jq dx<fi^ z (x, fify/x 
4.69, which is close to that of the original CZ-model 



In Ref. [54|, it was shown that by taking the asymptotic DA and considering the suppres- 
sion from the NLO contribution, the value of Q 2 F n ^(Q 2 ) at Q 2 = 8.0 GeV 2 is somewhat 
smaller than the experimental value (16.7±2.5±0.4)-10~ 2 GeV, i.e. Q 2 F 7T1 (Q 2 )\q2 =8 GeV 2 ~ 
0.115 GeV. However it should be pointed out that Ref. [54j| is only a simplified analysis, since 
the effects caused by the /cr-dependence in both the pion wavefunction and the hard scatter- 
ing amplitude, and the contribution from F^ V \Q 2 ) have not been taking into consideration. 
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As shown in Fig. (|7J) , even though Ff£j v \Q 2 ) is suppressed to Ff^\Q 2 ) in high Q 2 region, it 
has sizable contribution in the intermediate Q 2 region, e.g. it is about 12% at Q 2 = 8.0 GeV 2 
for the case of m — 0.30 GeV. 



D. A simple discussion of the NLO correction 

Up to NLO, by taking the square of the factorization scale y? F to be Q 2 , the pion-photon 
transition form factor F^^Q 2 ) can be schematically written as 

FAQ 2 ) = FiViQ 2 ) (l - 6 ■ + F^ V \Q% (40) 

where 5 is a parameter that is to be determined by the behavior of the pion wavefunction/DA 
and the detail form of the hard scattering amplitude. It is a natural choice to take fi 2 F = Q 2 , 
which directly eliminates the /^-dependence in the terms that determine how much of the 



collinear term is absorbed into the distribution amplitude lol.ll5(| . In the literature, the value 
of 5 has been estimated to be: 8 — | for asymptotic DA {9, 10, 18] and 5 — ^ for DA in CZ 



form [9;]. It should be noted that to be consistent with our present calculation technology, 
we need a full NLO result, in which all the effects caused by the transverse-momentum 
dependence in the hard-scattering amplitude and the wavefunction and by the Sudakov 
factor should be fully taken into consideration. However such a full NLO calculation is not 
available at the present. The value of 5 will be decreased by considering the transverse- 
momentum dependence in the hard-scattering amplitude and the wavefunction (a naive 



11|). Since our model wavefunction is close 



discussion for this point can be found in Refs. 9 
to asymptotic-like one, we simply take 5 = | to do our discussion. 

Under the condition of 8 = |, it can be found that the best fit of Q 2 F n ^(Q 2 ) for the case 
of the BHL-like wavefunction (j^3"j) is at m ~ 0.32 GeV. If taking m = 0.32 ± 0.10 GeV, we 
obtain ^(/W 2 ,) = — 0.02+o;gg, where fi ~ 2 GeV. Furthermore, if taking a 2 (l GeV 2 ) > as 
an extra constraint, we find that m must be in the region of [0.23, 0.30] GeV. Or inversely, 
we have a 2 (l GeV 2 ) G [0,0.06] for m 6 [0.23 GeV, 0.30 GeV]. 



V. SUMMARY 



In this paper, we have given a careful analysis of the pion-photon transition form factor 
F-n-yiQ 2 ) involving the transverse momentum corrections with the present CLEO experimen- 
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tal data, in which the contributions beyond the leading Fock state have been taken into 
consideration. As is well-known, the leading Fock-state contribution dominates the pion- 
photon transition from factor F nl (Q 2 ) for large Q 2 region and it gives only half contribution 
to F ny (0) as one extends it to Q 2 = 0. One should consider the higher Fock states' contri- 
bution to F^(0) at the present experimental Q 2 region. We have constructed a phenomeno- 
logical expression to estimate the contributions beyond the leading Fock state based on the 
limiting behavior of F^ V \Q 2 ) at Q 2 — > 0. The calculated results favor the asymptotic-like 
behavior by comparing the predictions from different type of the model wavefunctions with 
the experimental data. 

On the other hand, the present CLEO data provides the important information of the 
pion DA as one has a complete expression for the pion-photon transition form factor relates 
two photons with one pion. Our expression for F ni (Q 2 ) only involves a single pion DA. Thus, 
comparing the calculated results of F^(Q 2 ) by taking the BHL-like pion wavefunction with 
the CLEO data one can extract some useful information of the pionic leading twist-2 DA. 
Our analysis shows that (1) the probability of finding the leading Fock state in the pion is 
less than one, i.e. P q q = 56% and (r 2 +) qq = (0.33 fm) 2 with m = 0.30 GeV. This means that 
the leading Fock state is more compact in the pion and it is necessary to take the higher 
Fock states into account to give full estimation of the pion-photon transition form factor 
and other exclusive processes. (2) under the region of m e [0.20 GeV, 0.40 GeV], we have 
the DA moments: a 2 {nl) = 0.002±^054, a^o) = -0.02218;^ and all of higher moments. 
Such result is helpful to understand other exclusive processes involving the pion. 
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